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O . Abstract 

^ . 

^ , We study the chiral-odd transversity spin-dependent quark distribution function 



hi{x) of the nucleon in a constituent quark model. The twist-2 structure functions, 
O , fi{x), gi{x) and hi{x) are calculated within the diquark spectator approximation. 

Whereas an inequality fi{x) > hi{x) > gi{x) holds with the interaction between 
^ , quark and diquark being scalar, the axial-vector effective quark-diquark interaction, 

' which contributes to the ci-quark distribution, does not lead to such a simple relation. 

We find that hi{x) for the (i-quark becomes somewhat smaller than gf{x), when we fix 
the model parameter to reproduce known other structure functions. We also include 
corrections due to the non-trivial structure of the constituent quark, which is mod- 
eled by the Goldstone boson dressing. This improves agreements of fi{x) and gi{x) 
with experiments, and brings further reduction of hf{x) distribution. Consequences 
for semi-inclusive experiments are also discussed. 
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1 Introduction 

Study of the nucleon spin structure triggered by the unexpected EMC resuh 
becomes one of the most challenging issues in the hadron physics. The spin structure 
function obtained by the polarized lepton-hadron deep inelastic scattering tells us a 
lack of our understandings of the quark and gluon dynamics in the nucleon. iFiom 
the lepton-hadron deep inelastic scattering, one can obtain the longitudinal spin 
structure function (71(0;) and the transverse spin structure function, g2{x), which 
contains information on twist-3 pieces]^]. These structure functions reveal the helicity 
distribution function as well as the quark-gluon correlation in the nucleon|^. Further 
approach using the semi-inclusive process is now planed, which makes it possible to 
know detailed decomposition of the nucleon spin structure. 

Recently, chiral-odd structure function is discussed as complementary tool to 
study the structure of the nucleon ^. In the lepton-hadron scattering, in which the 
chirality is conserved, only the chiral-even structure functions can be observed. On 
the other hand, in the Drell-Yan process, the chirality changing process is possible, 
and thus the chiral-odd structure function is expected to be measured. 

In terms of the quark-quark density matrix, the twist-2 structure functions of the 
nucleon are expressed as[^], 
r fix 

/ -e''^ {P, S\miM^n)\P, S) = 2 + h.t.] (1) 

J — e'^^ (P, S\mi,l'H^n)\P, S) = 2 [g,{x)P,{S ■ n) + h.t.] (2) 

J —e'^- (P, 5|^(0)a^,^7'^(Ari)|P, S) = 2 [h{x){S^,P, - S^^P^) + h.t.] (3) 
(h.t. denotes the higher twist contributions.) 

Here, we followed the definition of Ref. 0. is the usual spin-independent quark 

distribution, and gi{x) the quark helicity distribution function in the longitudinally 
polarized nucleon. These two distribution functions are chiral-even, and can be mea- 
sured by the lepton-hadron scattering. The transversity spin distribution function 
hi{x) is chiral-odd, in which we are mainly interested in this paper. The hi{x) struc- 
ture function corresponds to a target helicity-flip amplitude in the helicity basis, and 
thus does not have simple partonic probabilistic interpretation. As pointed out in 
Ref. 0, however, the transversity spin distribution can be understood as difference 
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of numbers of valence quarks with eigenvalues +1 and —1 of the transverse Pauli- 
Lubanski operator in the transversely polarized nucleon. 

Our aim of this paper is to present hi{x) using low energy effective quark model, 
which is already used to calculate known structure functions fi{x) and gi{x). In 
particular, we discuss the flavor dependence of the transversity spin distribution 
hi{x) in detail, compared with other twist-2 quark distributions. We evaluate quark 
distribution functions fi{x), gi{x) and hi{x) of the nucleon based on a constituent 
quark picture, which is successful to describe the low energy hadron properties. Here, 
current quarks are assumed to acquire their dynamical masses due to spontaneous 
breakdown of the chiral symmetry following to the Nambu and Jona-Lasinio model[0. 



We take the model developed in Ref. P, |10| to calculate the structure functions, 
where the nucleon is described as a bound state of the quark and diquark. 

For effective interaction between quark and diquark, we take both scalar and 
axial-vector types from the symmetry consideration. In principle, we could solve a 
bound state equation of the quark-diquark for the nucleon. However, it needs difficult 



procedure, and is first recently done by Kusaka et al. for only the scalar channel [[IT |. 



Here, we simply assume reasonable forms of the quark-diquark vertex functions for 
the scalar and axial-vector channels so as to reproduce the observed structure function 
fi{x) and (71 (x), and try to describe hi{x) structure function. 

There are several works to investigate the transversity spin distribution function 
hi{x) within the quark models]^, |T2|, |TB|, |T^. Since hi{x) is defined by the chirality- 



violating operator, the transversity spin structure function may depend on the chiral- 



odd operators of the theory, e.g. the quark condensate of the QCD vacuum[|T2|. 

Recently, Soffer has proposed that such twist-2 quark distribution functions, /i (x) , 
gi{x) and are subject to Soffer' s inequality, which is derived from a general 

argument [1T5|, |16|. 



fi{x)+g^{x)>2\h^{x)\ (4) 

We discuss behavior of this inequality in our model. 

In particular, the simple relativistic quark model such as the bag model0, where 
the relativistic quarks independently move inside the confinement region, leads to a 
relation; 

/i(x) + (7i(x) = 2|/ii(x)| (5) 
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Mx)\ > Mx)\ (6) 

The first equation (^) means a saturation of tfie Soffer's inequality in tlie model, 
discussed by Goldstein et al. JlBl . The second inequality comes from roles of the lower 



j[dr] 
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j[dr] 




3^. 



component of the Dirac spinor. Indeed, 1st moments of gi{x) and hi{x), axial charge 
and tensor charge, are given by 

Aq = 

6q - 

where F and G are upper and lower radial wave functions of the Dirac spinor. Eq. 
is a universal relation of the relativistic quark potential model independent of the 
model parameters. 

However, recent lattice QCD simulation indicates the 1-st moment of the |/ii(a;)| 
is larger than that of |5'i(a;)| for the rf-quark distribution function ||17||. In addition. 



the QCD sum rule study suggests a very small value of the (i-quark tensor charge P|. 

In the present study, we take both scalar and axial-vector quark-diquark-nucleon 
vertices, and examine such relations. When the nucleon is a bound state of the 
quark and the scalar diquark, we exactly obtain + gi{x) = 2hi{x) and hi{x) > 
5'i(^)[|Hl) "which is quite similar with the bag model calculations. However, for the 
axial-vector vertex case, which determines the (i-quark distribution function, such 
simple relations are not maintained. We find \hf{x)\ < \gf{x) \ or they are comparable 
at least, if we fix the model parameters to reproduce the other structure functions. 

We also study corrections to the quark distributions from dressed structure of 
the constituent quark. Constituent quarks are assumed to be quasi-particles of the 
QCD vacuum, and have non-trivial structure. We introduce the Goldstone boson 



cloud, 7i,K,r], around the constituent quark studied in Ref. [|T9|. Such a dressing 
produces crucial effects on the structure function, namely, reduces a probability to 
find a bare quark state and also changes the spin structure by emitting the GS 
boson to relative P-wave state. Combining these contributions with the results of 
quark-diquark model, we show that shape of the transversity spin distributions are 
considerably different from those of the helicity distributions. 

This paper is organized as follows. In Section 2, we calculate the chiral-odd type 
forward scattering amplitude introduced in Ref. , which provides the transversity 
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spin distribution hi{x). Relations of three structure function are examined in some 
detail. In Section 3, numerical results are given performing the evolution of the 
distribution functions. and gi{x) are compared with the experimental data, and 

hi{x) is presented for forthcoming experiment. We introduce the dressed structure 
of the constituent structure due to the Goldstone boson cloud, developed recently by 
Suzuki and Weise|T^. Inclusion of the GS boson dressing improves the agreements 
substantially. Final section is devoted to summary and discussions. 



2 Calculations of the twist-2 quark distribution functions 

Purpose of the present paper is to study the structure function measured at high 
energies by means of the low energy quark model of the nucleon. Basic procedure 
to connect calculations of the quark model with the high energy experimental data 
is summarized as follows We evaluate twist-2 matrix elements of the structure 
function within the effective quark model at the scale fi ~ IGeV, where the effective 
models are supposed to work. Calculated structure function has no physical meaning 
at this scale /x, but plays a role of the boundary condition for the QCD evolution 
equation. With the help of the perturbative QCD, the results are evolved from the 
low energy model scale /z to the high momentum scale, at which the experimental 
data exist. Higher twist contributions can be neglected, since we are interested in the 
Bjorken limit oo. Comparison of the model calculations with the experimental 

data reveals non-perturbative dynamics in the deep inelastic experiments. 

Let us start calculating nucleon structure functions, fi{x), gi{x) and hi{x). We 
adopt the phenomenological quark-diquark model for the nucleon p, |^, In this 
approach, the virtual photon-nucleon forward scattering amplitude is illustrated in 
Fig.l. Here, the constituent quark is struck out by the high momentum virtual photon 
with a residual diquark being spectator. We treat diquarks as only the spectators, 
though the structure of the diquark itself might be crucial when we go beyond the 
diquark spectator approximation Regarding the nucleon-quark-diquark effective 
vertex, we deal with both scalar and axial-vector types, which are suggested by the 
success of the SU (6) spin-flavor symmetry approach. 

Vs = IMp') (7) 
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Figure 1: The forward scattering amplitude in the quark-diquark model 

The thick and thin solid lines denote the nucleon and the quark, respectively. The 
spectator diquark is depicted by the dashed line. The circle and box indicate the 
external currents, which are the axial- vector and scalar currents for hi{x). For fi{x) 
and gi{x) cases the current probes the nucleon structure as the vector and axial- 
vector, respectively. 

v{f = r^'Mp') (8) 

where (psip"^) and (pvip'^) are momentum cutoff functions, specified later. 

In the following studies, we use standard definitions for the quark distribution 
functions @, 

Fiix) = lel [mx) + ffix)] 

where Fi{x) is the standard electro-production structure function, and Ca is the elec- 
tric charge of the quark. Summation a is over quarks for all flavors. 

At first, we show basic steps to derive the spin independent distribution 
for the scalar channel. Explicit expressions for fi{x) and gi{x) are already given in 



literature p, 10 1. We define the virtual photon- nucleon forward scattering amplitude 



T,,{q,P,s) = zJ -0-^e'i<{P,s\TJ,{OJM\P,s) (9) 

where is the vector current. Here, q is the momentum of the virtual photon, and p 
the nucleon momentum with the mass = M^. s is the nucleon spin, which satisfies 
p ■ s = and = —1. We also introduce variables, v = P ■ q/M and the Bjorken-x 
X = — g^/2Mz/ = Q"^ /2Mv. By virtue of the optical theorem, this forward scattering 
amplitude is related with the hadronic tensor, which provides the structure function. 
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-Tr 



V^{k)S{k)^^S{k + q)^,S{k)V,{k)D,{k2){jP + M){l+^'' /,) (10) 



In this model, the matrix element is calculated as, 
d^p 1, 
(2^2 

where momenta delivered by the quark and diquark are k and k2 = P—k, respectively. 
Here, the quark propagator S{p) with the constituent quark mass m is written as, 

S{p) = (11) 
p — m 

The scalar (spin-0) diquark propagator is given by, 
1 



Ds{p) 



9 2 ' 

— TTLj-) 



and for the axial-vector (spin-1) diquark 
1 



Dv{p) 



(12) 



(13) 



— m 

Here, m^) is masses of the diquarks 
Using the identity, 

one can get the spin independent and dependent distribution functions after suitable 
projection^]. For instance, using the projector Vap = \ -^PaPp — Qap , we can get 
the unpolarized quark distribution /i(x). 

We introduce the light-cone variable, = x°±x^. Non- vanishing contribution to 
the imaginary part of the forward scattering amplitude arises from poles of the struck 
quark propagator [(/c + g)^ — m^] and the spectator diquark propagator [kl — mjj]. The 
former is reduced to a condition that the Bjorken-x is identified with the longitudinal 
momentum fraction of the parton. 



m 



6[ik + qy 
The latter one yields 



2P-q 



-6 



k+_ 



X 



5[kl - ml] = 5[{1 - x)P+k^ -kl-m] 

Inserting the above results into Eq. ([10|) , one finds 

„ . . [ dk^dk'dPkT 

fs[x) 



X 



2(27r)4 (P - m2)2 
2'n5{k+ - xP+) 1 
2P -q 2 
^a. de 0|(A;2) 



2tx5 



(1 - x)P- 
Tr[(^ + m)^(^ + m)(f +M)] 



2^2 



87r2 (F - 

X [x(M2 + 2mM - ml) + - (1 - x)k'^] 



(14) 
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where 

kL. = -T^< + xM^ ■ (15) 
1 — X 

It is possible to rewrite Eq. (|lj) as 

with the transverse momentum kjn given by 

fc2 = (1 - x)xM^ - xm|, - (1 - x)k'^ 

Now we come back to the chiral-odd distribution function hi{x), which can not 
be probed by the electro-magnetic current, as mentioned in Section 1. loffe and 



Khodjiamirian introduced the following scattering amplitude |]T2|, which allows us to 

evaluate hi{ 

T,{q,P,s) = tJ dV-«l(P,s|T{J^5(O^(0) + J(O^M5(0)}|P,s) (17) 

where we choose the transverse nucleon spin s = (0, 1, 0, 0). This forward scattering 
amplitude is related with the hi{x) structure function. 

Tf,{q, P, s) = (^s^ - ^i(^) 

+ i^^l' " ~^^^) ' ^^^^^^^ + £iJ,uXaq\Sa{q ■ S)l2{x) (18) 

Hence, one can obtain the transversity spin distribution function through the optical 
theorem; 

hi{x) = --lmhi{x) (19) 

TT 

Straightforward calculations yield the following expression for the scalar spin-0 di- 
quark spectator pro cess [p!^; 

hsix)= I ^tS^^ZZ^ {xM + mf (20) 



where we use {k ■ sY = k^ = k^/2. It is interesting to note that resulting hs{x) is 
proportional to masses of the nucleon M and the constituent quark m, which are 
the chiral symmetry breaking parameter of the model. If we assume a 'un-realistic' 
chiral symmetric world, M = m = 0, the hs{x) distribution function vanishes. This 
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is different from fi{x) and gi{x), which are non-zero even if all the particles are 
mass-less. 

Using ( ]TT| ) and (|I^) and summing spin-1 diquark helicities, we also get the hi{x) 
for the axial- vector diquark spectator process. 



ml 



(21) 



It is not easy to see how hf{x) behaves in the chiral mass-less limit. Even if we take 
m = M = 0, a term ~ k^/ml survives and hy{x) 7^ 0. 

Before constructing the nucleon structure function with the explicit spin-flavor 
wave function, we shall discuss relations of these quark distribution functions and 
the Soffer's inequality |15|]. We show in Fig. 2 fs{x), gs{x) and hs{x) calculated by 
the quark-scalar diquark model. In this case, the relation fs{x) > hs{x) > gs{x) 
holds, independent of the model parameters, and the Soffer's inequahty is saturated 
fs{x) + gs{x) = '^hs{x). This behavior is similar with the bag model. 

The spin-1 diquark spectator process shows rather non-trivial behavior. The quark 
distribution functions are shown in the middle of Fig. 2. The Soffer's inequality is 
satisfied, but can not be saturated. Moreover, relative magnitude of gv{,x) and hv{x) 
depends on the choice of the parameters. Here, we use the parameter set, which is 
suitable to reproduce other quantities such as the quark momentum fraction and gA- 
We can not find any simple relation for hy{x) and gv{x)- For completeness, we plot 
f (x) + g{x) —2\h{x) \ in a lower part of Fig. 2. In the spin-0 diquark spectator case, this 
quantity becomes exactly zero, which means a saturation of Soffer's inequality. Result 
of the spin-1 diquark case also satisfies a positive condition for f{x) + g{x) — 2\h{x)\. 

In terms of these distribution functions, we construct the nucleon structure func- 
tions with the spin and flavor structure. Since we have already included the spin 
structure of the diquarks to calculate the distribution functions and summed up over 
diquark helicities, we need only flavor structure of the nucleon wave function, which 
is given by 

\p> = -^S{ud)u+ ^A{ud)u ^A{uu)d, (22) 

v2 v6 v3 

where u,d denote the up-, and down-quarks, and S{ij) and A{ij) are the spin-0 
and spin-1 diquarks, respectively. Here we assume the S'f/(4) spin-flavor symmetry. 
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Figure 2: Quark distribution functions in the quark-diquark model of the nucleon 

The upper figure shows the results of the spin-0 diquark spectator process (scalar 
channel), and the middle one the results of the spin-1 diquark spectator (axial- vector 
channel). In each figure, the quark distributions f\{x), gi{x) and hi{x) are shown 
by the solid, dashed and dash-dotted curves, respectively. The lower figure shows 
fi{x) + gi{x) — 2\hi{x)\ for the scalar and axial- vector channel by the solid and 
dashed curves, respectively. 
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Quark distribution functions ^"(x) and qf{x) with q = f,g,h are expressed as 

q^,{x) = ^qs{x) + ^qy{x) , (23) 

qfix) = qv{x) . (24) 

Here qs and qv are corresponding quark distribution functions presented in Fig. 2 
with the spin-0 and spin-1 diquark being spectators, respectively. 

As easily seen in Eqs. (p3|j2^) , dominant contribution to the w-quark distribution 
comes from the spin-0 diquark spectator process, and the ci-quark distribution is 
governed by the spin-1 diquark process. Recall that Eq. ( PPD implies the transversity 
distribution vanishes when m = M = 0, hi{x) distribution would become very small, 
if the chiral symmetry were not spontaneously broken. 

Now we fix the model parameters with physically motivated constraints. The spin- 
diquark is expected to be tightly correlated quark-quark state due to the color spin- 
dependent interaction, as discussed in several works[p, |^. In particular, experimental 
data of the non-leptonic hyperon week decay require the strong correlation in the spin- 
diquark, which resolves the long standing problem of the AJ =1/2 rule[^, 

Here, we take into account differences between the spin-0 and spin-1 diquarks in 
the following ways. We assume that the mass difference of the diquarks is a few hun- 
dreds MeV from a consideration of the N-A mass splitting. Effective size of the spin-1 
diquark is expected to be larger than that of the spin-0 diquark, because it is assumed 
to be weakly bound state of quarks. We introduce the following parametrization for 
the quark-diquark vertex functions [|To|] . 



^sik) = 9SJ1-2 (25) 



(A| - k^y 

m? — k"^ 

Here coupling constants gsiQv are determined by the normalization of the unpolar- 
ized distribution fi{x). These functions ensure the correct behavior of the quark 
distributions at large x, /"(x) ~ (1 — x)^ and ff{x) ~ (1 — x)^, which are consistent 



with the standard parametrization ||25|. On the other hand, such large-x behavior does 



not agree with the parametrization of Brodsky et al.^6\. Recent reanalysis of the 



experimental data by Melnitchouk and Thomas ||27||, who take into account the Fermi 
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motion, nuclear binding and nucleon off-shell effects in the deuteron, is consistent 
with the QCD inspired fit by Brodsky et al.ER. 

Taking into account above conditions, we determine the model parameters by the 
experimental data. We take masses of the constituents, m = 450MeV, ms = 700MeV, 
and ruv = SOOMeV for the quark and diquarks. Values of the cutoff for the quark- 
diquark vertices are A5 = 0.49GeV and Ay = 0.66GeV. 

The resulting momentum fraction carried by the u and d quarks is found to be 
0.57 and 0.24 respectively, with the remaining 20% being left for the gluons and sea 
quarks. The spin fractions of the nucleon become Am = 0.95 Ad = —0.30, which 
yield the nucleon axial- vector coupling qa = 1.25. Calculated total nucleon spin is 
S = 0.65, which is similar with relativistic quark model calculations and larger than 
the empirical value. We also get the tensor charge, 6u = 1.17 and 6d = —0.26. 
We have a relation Au < 6u independent of the choice of the parameters. The d- 
quark case is more subtle as mentioned above. We get \6d\ < \Ad\, which seems 
to be consistent with the recent lattice result and disagrees with the bag model 
calculation [||. 

Such a tendency is quite different from the bag model result, which gives the 
universal inequality 6q > Aq (j^). In the relativistic quark potential model such as 
the bag model, quarks move inside the potential independently, and hence there is 
no correlation between the spin of the struck quark and spin structure of remain- 
ing spectator 2-quark, except for the trivial spin-isospin factor of the nucleon wave 
function. On the other hand, in a bound state of the quark and spin-1 diquark of 
the nucleon, helicities of the struck quark and spectator spin-1 diquark are correlated 
depending on the dynamics of the vertex function. Hence, it is rather reasonable that 
magnitudes of hf{x) and gf{x) are parameter dependent. 



3 Momentum evolution and comparison with experiments 

We have obtained the quark distribution functions in the nucleon at the low energy 
model scale fi. In order to compare then with experimental data, we carry out the 

evolution of the distribution function with help of the perturbative QCD. It can 
be done by using the Altarelli-Parisi equation for fi{x) and (71 (a;) . As for the 
transversity distribution, the Altarelli-Parisi splitting kernel is developed by Artru 
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and Mekhfip 



dhi{x,Q'^) OisiQ"^) f dy „ ix 



d{\nQ^) 277 I y ' \y 

The splitting function Pt is given by. 



Pt\-] h,(y,Q') (27) 



2 + U{z-l] 



f28) 



(1-2;)+ 2 

which differs from the fi{x) and gi{x) cases. Due to the current conservation, the 
1st moments of fi{x) and gi{x) are unchanged by the QCD evolution. However, the 
1st moment of hi{x), tensor charge, decreases with increasing the momentum scale. 

We take the low energy model scale /i^ = 0.2GeV^ with the QCD cutoff A = 
0.25GeV, which are adopted in Ref. |2^, and use the leading order (LO) evolution 
equations. As discussed in section 2, the low energy effective quark model should 
provide quark distributions at a typical hadronic scale ~ IGeV. However, if we 
adopted /i = IGeV as the low energy scale, the resulting distribution function would 
disagree with the experiment. In order to reproduce the data, the low energy scale 
must be smaller than IGeV, typically ~ O.SGeV. It is difficult to understand why 
such a discrepancy arises at the moment. Thus, we treat the starting scale /i as a free 
parameter so as to reproduce the experiments. Inclusion of the non-trivial structure 
of the constituent quarks might be crucial to study the origin of this discrepancy 0], 
which will be discussed below. 

Use of the perturbative QCD below IGeV seems to be questionable. However, 
inclusion of next-to-leading order (NLO) corrections modifies about 10 ~ 20% of the 
LO evolution result, which will be briefly discussed later. Of course, if we take much 
lower scale, e.g. ~ O.lGeV^, then difference of the LO and NLO results becomes very 
large, and the perturbative evolution is no longer reliable. All the results presented 
below are evolved to = lOGeV^. When we evolve the tensor charge from /i^ = 
0.2GeV^ to a few GeV^ scale, the reduction of the tensor charge is less than 20%. 

We first show /"(x) and ff{x) with the CTEQ experimental parametrization|25 
in Fig.3 at = lOGeV^ by thick-dashed and thin-dashed curves. The resulting dis- 
tribution at X ~ 0.3 is about 40% larger than the experimental fit. In the framework 
of the quark-diquark model (and most effective models), it is quite difficult to get 
better agreements within the reasonable parameters of the model. In the next sec- 
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- - xfi"(x) 

- - Xfid(x) 

— xfi"(x) with CQ 

— xfid(x) with CQ 
o CTEQ(u) 

o CTEQ(d) 




Figure 3: The unpolarized valence quark distribution of the nucleon 

Unpolarized valence u and d quark distributions at = lOGeV^. Thick-dashed and 
thin-dashed curves indicate results with bare constituent quarks for u and d distribu- 
tions, respectively. Results with the GS boson dressing are shown by thick-solid and 
thin-solid curves for u and d quarks. The circle denotes the CTEQ4 parametrization 
for the valence w-quark distribution, and the boxes for the ci-quark| 



tion, we will introduce dressing of the constituent quark due to the Goldstone boson 
cloud, which greatly improves the numerical results. 



Here, we comment on the work by Kulagin et al. [|T0[. Procedure to obtain the 



unpolarized structure function presented here is the same as one of Ref. |T0 |. 

They also estimate contribution from structure of the constituent quark itself, which 
is described by the pion dressing at the middle-x region and the Regge exchange 
contribution at the small-x. Combining those important contributions with the bare 
quark-diquark model, they have obtained improved quark distributions as the input 
of the evolution at typical hadronic scale ~ IGeV, at which the use of the per- 
turbative QCD evolution is reasonably acceptable. They have reproduced observed 
structure function of the nucleon very well. 

We have analyzed our results for the spin independent as well as the spin depen- 



dent structure function with the same parameter set used in Ref. Unfortunately, 
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we find resulting spin dependent structure functions become too small compared with 
the experimental data. This is because they adopt large diquark masses niD > IGeV. 
To reproduce the nucleon axial- vector coupling g^, it is necessary to take smaller val- 
ues of diquark masses around 700MeV within this quark-diquark model, and thus 
impressive agreements of with the data are lost. Calculated distributions with 
our parameters also require the evolution from /i^ ~ 0.2GeV^ to get reasonable agree- 
ments with the data. 

The longitudinal spin structure function of the proton Gi{x) = | ffi'i (a^) + ^gfi^) 
and the neutron = | ^diix) + ^gf{x) are shown in Fig. 4. Our calculation 

reproduces the nucleon axial-vector coupling = 1.25, but the quark spin frac- 
tion Am = 0.95 and Ad = —0.30 are somewhat different from the empirical values 



Au = 0.82 ± 0.03 and Ad = —0.43 ± 0.03 [^. As a result, calculation overestimates 
gi{x) and the total nucleon spin fraction S = 0.65 compared with the empirical 
one S = 0.31 ± 0.07. Note that gf{x) becomes positive at the large x in this model, 
which is consistent with the QCD inspired parametrization of the helicity distribution 
functions by Brodsky et a/. p6[|. 

Before discussing the chiral-odd structure function, we point out the importance 
of ratio of the structure functions in order to extract the non-perturbative aspects of 
QCD from the deep inelastic data. As mentioned above, shapes of calculated quark 
distribution depend so strongly on the choice of the low energy scale that it is not 
easy to draw definite conclusions from the comparison of the model calculations with 
experiments. However, if we take a ratio of several distribution functions, the calcu- 
lated ratio itself is rather insensitive to the choice of the scale /i^, since ambiguities 
from the evolution are considerably canceled out by taking a ratio. Therefore, we can 
know how non-perturbative structure of the nucleon affects the distribution functions 
in the deep inelastic scattering. Also, by taking a ratio, we can check internal con- 
sistency of the model calculations. Results for fi{x)/ fi{x) and Ai = G\{x) / Fi{x) 
presented in Fig. 5 show remarkable agreements with the data. These agreements 
are consequence of the spin-flavor structure of the nucleon, which is incorporated by 
the quark-diquark model|§]. If we took un-physical values for diquark masses, e.g. 
rris > my, agreement would be lost. Correct spin- flavor correlation of the nucleon 
wave function leads to a good description of these ratios. 

We present and hf{x) in Fig. 6 at = lOGeV^. The -u-quark distribu- 
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Figure 4: The longitudinal spin structure function of the proton and the neutron 

Upper figure shows the proton spin structure function Gi{x), and the lower one 
the neutron G"(a;) at = lOGeV^. In each figure, result with the bare quark is 
shown by the dashed curves, and one with the GS boson cloud by the solid curves. 
Experimental data is taken from Ref. [Q. 
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Figure 5: Ratios of the structure functions with the experimental data 

fi{x) / fi{x) and Gf(x)/Ff(x) are shown in upper and lower figures, respectively. In 
each figure, calculated result is shown by the solid curve with the data||2^, 
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Figure 6: The transversity spin structure function hi{x) for the u and d quarks 

h'lix) and hfi x) are displayed by the thick sohd and dashed curves, respectively. For 
comparison, gi{x) and gf{x) are shown by the dotted and dash-dotted curves. 

tion shows similar behavior with fi{x) and gi{x), and is satisfied with /"(a;) > 
h^{x) > gi{x). On the other hand, shape of the hf{x) distribution differs from 
that of gf{x), and their magnitudes are similar. We will discuss the flavor depen- 
dence of the transversity spin distribution function in detail after introducing the GS 
boson dressing in the next section. 

We remark that the helicity distribution functions gi [x) becomes much larger than 
the transversity distributions hi{x) at the small-x region, x < 0.1. Remember that, 
at the low energy scale, is smaller than h^{x) for all x, as shown in Fig. 2. This 

smallness of hi{x, lOGeV^) at a; < 0.1 is due to the difference of the QCD evolution. 
Regarding the 1st moment of hi{x) structure function, the splitting kernel Eq. ( |28D 
for hi{x) does not give a significant difference, i.e. the reduction of the 1st moment is 
less than 20% by carrying out the evolution from fj? to lOGeV^. At small-x region, 
however, the evolution procedure generates a quite large difference between gi [x) 
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Figure 7: Goldstone boson corrections to the quark distribution function 



Fig.|^(a) indicates the Goldstone boson spectator process, and Fig.|^(b) probes the 
internal structure of the Goldstone boson with the constituent quark being spectator. 
The constituent quark and the Goldstone boson are depicted by the solid and dashed 
lines, respectively. 



and hi{x). As discussed in Ref. [^, such a difference produce great influence on the 



pp Drell-Yan experiment, which is planed at RHIC-SPIN. Transverse Drell-Yan spin 
asymmetry becomes substantially smaller than the longitudinal one An, which 
makes it very hard to measure the transversity spin distribution function hi{x) at 
RHIC. 



4 Corrections from the constituent quark structure 

In the previous section, we treat the constituent quark as a point-like particle. 
Here, we shall consider contributions from the Goldstone boson dressing of the con- 



stituent quark studied in Ref. ||19|. In the chiral quark model of Manohar and 
Georgi|3^, the constituent quark interacts with the Goldstone (GS) bosons, vr, K 
and rj, associated with the spontaneous chiral symmetry breaking. Within the one 
GS boson approximation, diagrams illustrated in Fig. 7 contribute to the structure 
function. The GS boson fluctuation around the constituent quark changes probabil- 
ity to find a bare quark state, and gives soft contribution to the quark distribution 
function. Also, it changes the CQ spin structure by emitting the GS boson into the 
P-wave state. 

The constituent quark-GS boson splitting functions giving corrections to the twist- 
2 structure functions are given by|[l9[|. 



piyh.,. = i, I i'4 .'j''rr^:'j (29) 



7 y{l-y) 
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j a/i 



1 {rrij — rriiyY — 



(30) 



[mi 



yO--y) 



m. 



2 - M2_^ 



(31) 



for fi{x), gi{x) and hi{x) distribution functions, respectively. Here, mi^mj.ma are 
masses of the j- constituent quarks and the pseudoscalar meson a, respectively. 



2 ^ "^j+^T 

i° y 



H is the invariant mass squared of the final state. We use the 

quark-GS boson coupling constant g = 3.76 as suggested by the effective theories^. 
These functions are calculated in the Infinite Momentum Frame, -f- exponential cutoff 
is used with Aqs = 1.9GeV, which is determined to reproduce the violation of the 
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g ^ g exp 



2 

mf 



M2 



2A|5 



(32) 



In terms of these splitting functions, the quark distribution qj{x) generated by 
the emission of the GS boson a from a parent quark distribution qi{x) (Fig. 7(a)) is 
expressed as; 

y " ' ' ' \y, 

The CQ spectator process in Fig. 7(b), which contributes to the small x region, is 
given by. 



Qji^) = I —Pja/iiy)qt 

J X V 



(33) 



qk{x) 



dyi dy2,. /a; \ / 

yi 1/2 \yij \y2 



(34) 



where Vg/a(x) is the quark distribution function in the GS bosons. The splitting 
function is obtained by a symmetric relation as Paj/i{x) = Pja/i{^ — x). Since we 
incorporate only the pseudoscalar mesons, this diagram never contributes to the spin- 
dependent processes. 

Combining the flavor structure of the GS bosons, one can rewrite the twist-2 
structure functions renormalized by the GS boson clouds as follows. 



f^ix)R = Zf^{x)+Pu.,d®ft + Vu/^®P.d/u®n + -Pu.lu®fl 



■'■It is possible to calculate the dressing in the covariant approach, but convolution breaking terms 
appear and the resulting formulae become much complicated. It is also possible to evaluate both 
bare distributions and dressing in the infinite momentum frame from the beginning. 
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+ ^K/. ® P. ® + ft) + Vu/K ® Pku/u ® fl 

+ \Pur,/u ® /l" + ^Vu/^ ®P^® Ul + /') (35) 
D do 

+^Ka ® ® (/r + ft) + ^K/, ® A ® (/r + ft) (36) 

Similar expressions can be written for the ci-quark case. For the spin- dependent 
distribution functions, 

g\{x)R = Zg'^ix) + Pu^/d 0gt + ^Pun/u ®gi+ \Pur,/u ® < (37) 

gi{x)R = Zgi{x) + Pd^i^ ®g'l + ^Pd./d (E) gf + ^Pd^/a ® gf (38) 

Here, ® expresses the convolution integral. Formulae for hi{x) distribution are ob- 
tained by just replacing g with h in Eqs. (pTl) , (|38D . The renormalization constant Z 
is given by, 

Z=l-l{P^)-{PK)-l{Pr,) 

where (Pa) is the 1st moment of the unpolarized splitting function. Within our model, 
Z is found to be 0.76. We emphasize that AP is somewhat negative or nearly zero. 



while 6P is positive |19|. Calculated first moments of the splitting functions for pion 
case are (P^r) = 0.13, (AP^) = —0.06 and {6Pt,) = 0.07. This sign difference of the 
helicity and transversity distributions causes considerable effects on the proceeding 
calculations. 

The results for the unpolarized valence quark distribution defined by fi{x) — fi{x) 
are presented in Fig. 3 by using Eqs. (|35| , p^). Here, we show the original result by the 
dashed curves, and the results with the CQ structure by the solid curves. The QCD 
evolution procedure is the same as before with /i^ = 0.2GeV^. The renormalization 
of the bare quark distribution brings down the peak of the distribution function for 
both u and d quarks, and the agreement around a; ~ 0.1 is considerably improved 
due to the GS boson dressing. 

The GS boson contribution to g^ and gf are different, as shown in Fig. 8. As 
a result of the renormalization, gi{x) distribution function decreases substantially, 
while gf{x) is almost unchanged. This fact is easily understood as follows. The 
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Figure 8: Helicity distribution functions gi{x) and gf{x) with the GS boson correc- 
tions 

Results with the dressing are shown by the thick sohd and dashed curves for u and 
d quarks, and ones without the dressing by thin curves. 

renormahzation factor Z always reduce absolute values of the distribution function, 
i. e. positive contribution to original (negative) gf{x). The largest correction to gf{x)ji, 
the second term of Eq. (|38D AP^n/u ® fl'i , gives a negative contribution, since AP^n/u 
is negative. Thus, cancelation of the renormahzation effect and the correction term 
keeps gf{x)R almost unchanged. Resulting spin fractions are tabulated in Table 1 
with the original values and lattice data. 

This is not the case for the transversity spin distribution. Results are shown in 
Fig. 9 with the original ones and gf{x)R for comparison. Both h^{x)fi and hf{x)fi 
decrease, particularly, the (i-quark transversity spin distribution becomes about a 
half of the result without the GS boson dressing. This is because the GS-boson 
splitting function 6P is positive for the transversity spin case, opposite to the helicity 
distribution case. Hence, for the (i-quark, the renormahzation as well as a correction 
term 5Pd-K/u®h^ are positive, which drive further reduction of the (i-quark transversity 
spin distribution. At a; ~ 0.1, the transversity spin distribution becomes the one-third 
of the helicity distribution. Corrected tensor charges are given in Table 1. 
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Figure 9: Transversity spin distributions and h'^{x) with the GS bosons 

Notations are same as those in Fig.^ We also show gf{x) with the GS boson dressing 
by the dotted curve. 



22 



It is very interesting and important to consider experimental possibilities to mea- 
sure such a flavor dependence of the transversity spin structure function. The Drell- 
Yan process, which is planned in RHIC-SPIN or HERA-N is not suitable to determine 
the flavor decomposition of the structure function, since all possible combination of 
q-q pair contributes to the fj^~^fi~ cross section. On the other hand, semi-inclusive 
deep inelastic scattering may be a good candidate to observe the flavor dependence 
of the spin structure function. For example, according to the work of Kotzinian and 
Mulders we stress the semi- inclusive deep inelastic scattering on the transversely 
polarized nucleon. Let us consider the vr production process (momentum Pt^) on the 
proton target (momentum P) with unpolarized lepton beam (momentum I). They 
have found the following weighted transverse asymmetry of the cross section is related 
with the transversity spin distribution of the nucleon. 

Jd<f)iJd'P^^^-^sm{<f)i + <j)i){da+-da-) 
^' = Jd^Jd^P.,ida-^ + da-) ^''^ 

where scaling variables are given by x = Q^/2P-q, y = P-q/P-l, and the momentum 
fraction of a produced pion z = P ■ P-j^/P ■ q. P^i. is the transverse momentum 
perpendicular to q, and azmitual angles 0'^^ are defined with respect to the lepton 
plane. Asymmetry of this cross section is generated by the so called Collins effect in 
the fragmentation process[^]. Eq. ( PU] ) is expressed with the quark distribution and 
fragmentation functions as, 

where h1{x) and /{'(x) are the quark distribution functions, and D^{z) and H^^{z) 
unpolarized and transversely polarized quark fragmentation functions from a quark 
flavor a to the pion. It is well known that the vr"*" production process is dominated 
by the w-quark fragmentation, and 7t~ by the d-quark. Also, since pions are isospin 
symmetric objects, the quark fragmentation functions are the same for -u — > vr^ and 
d ^ . Therefore, if we take a ratio of asymmetries for tt^ and tt^ production 
processes under above assumptions, we can arrive at the following expression. 

^ILm = MMZlM (41) 

Af ^ ' h\{x) ft{x) ^ ' 

where unknown 'transversely polarized quark fragmentation function' is canceled out. 
fi{x)/ff{x) is well determined so far, thus we can extract a ratio of the transversity 
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spin distribution functions of u and d quarks by counting the numbers of produced 
TT"*" and TT^. Using a similar technique with the polarized lepton beams on the longi- 
tudinally polarized proton, it may be possible to obtain a ratio of helicity distribution 
functions []3l| . 



We show in Fig. 10 hf{x)/h'^{x) using our calculated transversity distributions, 
where we can see clear flavor dependence of the spin-dependent structure functions. 
We also show gf{x)/g^{x) in our model, and one obtained by the parametrization 
of Brodsky et a/. Due to the reduction of hf{x), hf{x)/h^{x) is much smaller 
than gf[x)/gi{x) around x ~ 0.1. Calculated gf{x)/g'^ (x) is consistent with the 
the parametrization of Brodsky et a/. ||26|. The bag model calculation gives similar 
behavior for hf{x)/hi{x) and gf{x)/gi{x). Note that, in the naive S'f/(4) symmetry 
limit, hf{x)/h'^{x) = gf{x)/gi{x) = —1/4. Alternative experimental processes to 
measure hi{x) are now under investigation. 



5 Discussions 

We have studied the twist-2 quark distribution functions fi{x), gi{x) and hi{x) 
of the nucleon in terms of the phenomenological quark-diquark model. We deal 
with the scalar and axial-vector effective nucleon-quark-diquark vertices from the 
symmetry consideration. Resulting distribution functions show a saturation of the 
Soffer's inequality fi{x) + gi{x) = 2hi{x) when we use the scalar type of the quark- 
diquark interaction, whereas the axial-vector case is not saturated. Moreover, the 
spin distribution function satisfies a universal inequality gi{x) < hi{x) for the scalar 
vertex case, which is similar with result of the relativistic potential model like the bag 
model. On the other hand, relative size of these distribution functions depends on the 
model parameters in the axial-vector vertex case. With the model parameters fixed 
to reproduce the known observables, the quark momentum fraction and the nucleon 
axial- vector coupling gA-, we have found that is smaller than |5'i(a;)| or they 

are comparable at least. Since the axial- vector case contributes to only the (i-quark 
distribution from the SU{Q) symmetry, this behavior is consistent with the recent 



lattice QCD calculations and the QCD sum rule result [0], where the (i-quark 
tensor charge is substantially suppressed compared with the axial-charge. 

This difference simply comes from the spin structure of the quark-diquark vertices. 
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Figure 10: Flavor dependence of the spin dependent structure functions 

Ratios hf{x) /hi{x) and gf{x)/gi{x) are shown by the sohd and dashed curves. The 
QCD inspired parameterization p6| for gf{x)/g^{x) is also shown by the dotted curve. 
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We have checked other types of the Dirac structure for the vertex, and found, if the 
vertex is independent of the spin of the spectator diquark, e.g. scalar, pseudo-scalar, 
derivative coupling, universal relations Eqs. (^,^) hold. However, for the vector or 
axial-vector interactions, such relations are not realized, depending on the dynamics 
of the quark-diquark system. 

We have also estimated corrections due to the non-trivial structure of the con- 
stituent quark itself. The constituent quark and Goldstone bosons are believed to be 
most important degrees of freedom below the chiral symmetry breaking scale ~ IGeV. 
Corrections from the GS boson fluctuation can not be negligible, because it produces 
a renormalization of the CQ state as well as it changes the CQ spin structure by 
emitting the GS boson into the P-wave state relative to CQ. The probability to find 
a bare constituent quark decreases to about 70%, and the quark distribution func- 
tion becomes soft. Although corrections to the spin dependent parts are parameter 
dependent, the GS boson fluctuation gives positive contributions to the transversity 
spin distribution hi{x), while negative for the helicity distribution gi{x). This feature 
causes further reduction of the (i-quark tensor charge as shown in Table 1. 

Recently, experimental methods to test the nucleon transversity spin distribution 
are discussed by several authors P, |3^. The transversely polarized pp Drell- 



Yan process is planed at RHIC, but it is suggested that the transverse double spin 
asymmetry becomes quite small in the accessible kinematic region at RHICp9|. This 
is due to the suppression of hi{x) at x < 0.1 shown in Fig. 6. This is a general result 
of the effective quark model calculation, if one accepts calculations of the effective 
theory as quark distributions at the low energy scale /i ~ O.SGeV. The QCD evolution 
procedure from /i^ to the experimental scale generates a large difference between hi (x) 
and gi{x) at x < 0.1, even if the effective quark model provides the same inputs for 
hi{x) and gi{x). 

In Ref. 1^^, Artru et al. have given a constraint for the transversity spin distri- 
bution in the nucleon from the analysis of the single spin asymmetry in the pion 
production in the transversely polarized proton-proton collision p + p }^ tt + X. 
Their best fit gives 
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which is quite different from the simple SU{6) values, 

h'lix) _ 2 hf{x) _ 1 
7fM~ 3' 7fixj~~3 • 

We show our numerical result in Fig. 11. hf{x)/ ff{x) and gfix) / ff{x) are shown 
by the solid and dashed curves, and constraint from (^) by crosses. Our results 
are well consistent with the constraint (|42|) from the single pion production. The 
QCD parametrization of gf{x) / ff{x) is also shown by the circles Note that 

calculated distribution functions never violate the Soffer's inequality in any cases. 
Rather singular behavior of = —1 at a; = 1 is necessary to reproduce the 

observed large spin asymmetry of the tt" production. If we used the simple SU{Q) 
value hf{x) / ff{x) = —1/3, the resulting spin asymmetry would become about a half 
of the experimental value [^]. It is of great interest to investigate further approach 
to measure the hi{x) structure function, in particular its flavor dependence. 

There are some difficulties to calculate the structure function within the effective 
quark model. Most serious one is use of the perturbative QCD below 1 GeV. To 
obtain better agreements with the experimental data, the evolution from too small 
scale /i^ = 0.2 ~ O.SGeV^ is needed. We show in Fig. 12 obtained by LO and 

NLO evolution to clarify size of higher order corrections. Here, calculations with LO 
and NLO evolution are displayed by the solid and dashed curves, respectively. We 
use the same value of Aqcd to calculate both LO and NLO results for simplicity, 
though the LO and NLO values should be in general different as found in various 



QCD fits of the structure function |P4||. It is found that difference is about 20% of the 
original result, which is not significant. 

Whatever difference between the NLO and LO calculations is about 20%, validity 
of the perturbative QCD at such a scale is still questionable. It has recently been 
pointed out by Kulagin et a/.[|^ that inclusion of the constituent quark structure 
improves this difficulty. They have shown that, if one takes into account the inter- 
nal structure of the constituent quark, resulting distribution function becomes much 
softer. Using such improved distributions as inputs of the evolution at IGeV^, 
where the use of perturbative QCD can be justified, they have obtained remarkable 
agreements with the data. Our present work is extension of their work to the spin 
structure function. Although we have found that their model parameters gives too 
small values for the spin dependent properties like qa, such efforts contribute to right 
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Figure 11: Ratios of spin distribution functions hf{x)/ ff{x) and gfix)/ ff{x) 

hf{x)/ fi{x) and gf{x)/ff{x) are shown by the sohd and dashed curves, respectively. 
Best fit parameterization for hf{x)/ fi{x) extracted from the single spin asymmetry 
of the pion production in the polarized pp collision is also shown by the crosses. 
The circles denote the parameterization of Brodsky et al. for gf{x)/ff{x) 
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Figure 12: Comparison of LO and NLO QCD evolution 

Results for unpolarized w-quark distribution functions with the NLO and LO 
evolution are shown by the solid and dashed curves. 



direction to overcome this difficulty. 

Our aim is to investigate the non-perturbative nucleon structure from studies of 
the quark distributions measured in the high energy experiments. Obtained sup- 
pression of the (i-quark tensor charge is a key to understand how the real nucleon 
structure differs from the simple bag model-like structure. Future experimental efforts 
will clarify spin-flavor structure of the nucleon. 
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Table Captions 

Table 1 

The first moments of the helicity and transversity distribution functions (at /x^). 
The axial charges Am and Ad are shown in the second and third columns, and corre- 
sponding tensor charges are given in the forth and fifth columns, respectively. Results 
of with the bare CQ and dressed CQ are displayed in the second and third raws. Ex- 
perimental data, lattice calculations and QCD sum rule results are shown in the third, 
forth and fifth raws, respectively. 

Table 1 





Am 


Ad 


6u 


6d 


Bare 

With CQ 
Exp |28| 

Lattice|17| 
Sum rule 14 


0.95 
0.73 
0.83 ±0.03 
0.76 


-0.30 
-0.25 
-0.43 ±0.03 
-0.35 


1.17 
0.92 

0.84 
1.33 ±0.53 


-0.26 
-0.15 

-0.23 
-0.04 ±0.02 
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